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Abstract 

We consider a mechanism of generation of huge waves by multi-soliton resonant inter- 
actions. A non-stationary wave amplification phenomenon is found in some exact solutions 
of the Kadomtsev-Petviashvili (KP) equation. The mechanism proposed here explains the 
character of extreme waves and of those in Tsunami. 

1 Introduction 

The phenomenon of particularly high and steep waves on the sea surface is one of the 
most dangerous events. The large waves such as tsunami and freak (rogue) wave have a 
significant impact on the safety of people and infrastructure and are responsible for the 
erosion of coastlines and sea bottoms and the biological environment. For example, the 
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Indian ocean tsunami of 26 December 2004 caused an estimated 250,000 or more deaths 
and extensive damage due to run-up, landward inundation, and wave-structure interactions. 
The prediction of extreme waves is an important task for human being. Understanding the 
physics of large extreme waves may save lives. 

Freak waves occur much more frequently than it might be expected from surface wave 
statistics whereas they are particularly steep {TJ. Freak waves may happen in both deep 
water and shallow water. Although several physical mechanisms of generation of freak 
waves in deep water have been proposed, the mechanism in shallow water is few Q. 

Recently, it has been established that a moving pressure disturbance can generate soli- 
tary waves also in open sea areas. This mechanism may occur in relatively shallow areas 
with heavy fast ferry traffic where soliton-like disturbances frequently occur. In a fatal acci- 
dent which occured in Harwich, a port on England's east coast, in July 1999, one surviving 
victim reported that the wave looked like 'the white cliffs of Dover'. Research carried out 
in Europe shows that the soliton produced by a fast ferry was probably responsible for the 
disaster. This extreme wave is so-called 'solitary killer' wave (3]|4l|5l. 

Solitary waves in shallow water have been studied since the 19th century. The theories 
of asymptotic expansion and integrable system have brought the remarkable progress to this 
research field. With the assumption of weak nonlinearity the dynamics of solitary waves 
in shallow water is well described by the Korteweg-de Vries KdV) equation. Moreover, 
for weakly two-dimensional cases, that is, for the cases that the scale of variation in the 
direction normal to the propagation direction is much longer than that in the propagation 
direction, the fundamental equations may be reduced to the Kadomtsev-Petviashvili (KP) 
equnation. Both of the equations are integrable and it makes the precise analysis of the 
wave motion possible. 
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Peterson et al. studied the amplitude of two soliton solutions of the KP equation. It was 
pointed out that the interaction of two solitary waves may be one of mechanism of genera- 
tion of extreme waves (6l|71[8l. In the case of resonant Y-shape soliton which was found by 
Miles, the maximum amplitude of a solitary wave can reach four times the amplitude of an 
incoming solitary wave @[lOl. The weak point of this theory is that the interaction pattern 
of 2-soliton is stationary and does not describe non-stationary phenomena. If a high wave 
is generated by soliton resonance once, the high wave can not disappear within finite time. 
There is a gap to describe the generation of extreme wave using 2-soliton solutions of the 
KP equationSee ifTTTl . 

Recently, Biondini and Kodama studied a class of exact solutions of the KP equation 
and found that the interaction patterns of solitons in the class have non-stationary web-like 
structures [12]. Such solutions are resonant-type soliton solutions which are the general- 
ization of the Y-shape resonant soliton solution which was found by Miles. These general 
resonant solutions are simply expressed by a Wronskian. Furthermore, Kodama gave a 
general formulation of the classification of interaction patterns of the N-soliton solutions of 
the KP equation d (See also ifUl ). 

In this paper, we propose a mechanism of generation of huge waves by multi-soliton 
resonant interactions. The key of analysis is the Wronskian solution of the KP equation. 

2 Physical derivation of the Kadomtsev-Petviashvili 
equation 

Let us consider an object region of the ocean as an inviscid liquid layer of uniform undis- 
turbed depth h while the influence of the open air is assumed to be negligibly small. Assume 
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the plane z* = of the Cartesian coordinates (x* ,y*,z*) coincides with the fiat bottom of 
the region and z* = h + (x* ,y* ,t*) denotes the free surface. Here, t* is time. Let us 
denote the velocity components along the axes x*,y*,z* by u*(x* ,y* ,z* ,t*), v*(x* ,y* ,z* ,t*) 
and w* (x* ,y* ,z* ,t*), respectively. Assume the fluid motion is irrotational, then the fluid ve- 
locity components are written as u* = <J>** , v* = <&*„ , and w* = <!>*» in terms of the velocity 
potential <t>*(x* ,y* ,z* ,t*). Here we consider the long waves of small but finite amplitude. 
Hence, we make the nondimensinalization as follows: 

(x*,y*) = £(x,y), z* = hz, f = -t, C = at;, <I>* = (1) 

c h 

where £ is the characteristic horizontal length, a the characteritic amplitude, g the gravita- 
tional accelaration, and cq = yfgh the linear long-wave phase speed. Then the governing 
equations and the boundary conditions are: 

j3A* + * a = 0, (0< Z <1 + «C), (2) 

<*> z = 0, (z = 0), (3) 

^ + |(VO) 2 + ^<I>2 + C = 0, (z=l + aC), (4) 

£ + aVS-VC-^ = 0, (z = l + aQ, (5) 

where the subscript variables denote the partial differentiations with respect to the variables 
and V = (d/dx,d/dy), A = d 2 /dx 2 + d 2 /dy 2 , a = a/h, and j8 = (h/£) 2 . Here we consider 
weakly nonlinear long waves where weak nonlinearity and weak dispersion balance out, 
thus we assume a<l,^<l and j8 = 0(a) (hereafter we assume e = a = j3, that is, e = 
a/h and £ = h/y/e). Then, from (f2]) and ([3]), the velocity potential &(x,y,z,t) is represented 
as follows: 

°° Jim 

®(x,y,z,t) = £ (-eA) m 4>{x,y,t)-±— , (6) 

m=0 v Zm / • 



where (j>(x,y,t) is the value of the velocity potential at the bottom. Substituting © into the 
boundary conditions (@]) and ([5]), we obtain the coupled equations for the functions and £ 

<k + C + ^V0.V0-i £ A0 f = O(£ 2 ), (7) 

C ; + A4> + £V-(CV0)-i £ A 2 = 0(e 2 ). (8) 
o 

For unidirectional propagation (in the positive x direction) of waves, if we introduce the 
new variables t, =x — t, a = et and assume the expansions 

C = Sb(€ I <T) + eCi(€ I <T) + ---, = 0o(£,<7) + e<M£, <*) + ••■, (9) 
we obtain the KdV equation for C$ : 

As is well known, this equation describes the dynamics of one dimensional soliton (line 
soliton). Miles formulated the interaction of two solitons propagaing in different directions 
ni and 112 by means of a perturbation method |9J- Miles showed the followings: If the 
parameter k: = |(1 — ni • 112) = sin 2 \ff (2y is the angle between ni and m) is not small, 
in the lowest order of approximation £ can be represented by the superposition of two 
solitons propagating in the directions ni and 112 which are governed by the respective KdV 
equations. But, if K = 0(e), the perturbation expansion breaks down. Miles called the 
former the weak interaction and the latter the strong interaction. 

In the case of K = 0(e), the directions of propagation of two solitons are almost the 
same. If we choose the x axis as the main propagation direction, the scale of variation of 
the wave field in the y direction would be f?(e -1 / 2 ). Thus, introducing the variables 

£=x-t, rj=e 1/2 y, o = et, (11) 
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and expanding £ and <p as 



C = Co(£,T?,<x) + eCi(£,77,<7) + ---, =0o(£,T7, <r) + e0i(£, r], o) + -~, (12) 



we obtain the KP equation (HI [HI 



or 



V ^ + ?0 W + 3 W J + W ~ ( } 



In terms of the dimensional variables, this is written as 



d fdQ dQ 3c dQ 1 2 d 3 C *V c d 2 C * n « 



where £q = a£o- 

Here, if we introduce the following dimensionless variables: 



then (1151) becomes 

(-4itr + Mjfxx + 6wMjf)x + 3MFF =0. (17) 
Further, introducing T-function through 

u = ^L =2(logz) xx , (18) 
2h 

we obtain the following bilinear form from (fTTT ): 

(D x (-4D r +D|) + 3D F )T-T = 0, (19) 
where Dx, Dy and are the Hirota bilinear operators defined, for example, by 

D n x D™f.g=(d x -ds) n (d y -dy) m f(xj)g(y,y%= x ^ y . 
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Using the Hirota bilinear method ifTTl . the T-f unction of 1 -soli ton solution is given by 

T = l + exp(2@), © =K-X + QT + 0° = KX + LY + QT + 0° , (20) 

n=K { K2+3 S)> (21) 

and then the dimensionless elevation £ = Q/h is given by £ = ^(logr)xx = ^|-sech 2 ®. 
The T-function of 2-soliton solution is given by 

T = 1 +exp(2@i) +exp(2© 2 ) +Aexp(2@! +2© 2 ) , (22) 
©, = Ki • X + n t T + Qf = KiX + LjY + Q{T + &f , (23) 

n i = Ki (Kf + ^Y (i = l,2), (24) 

= 4( j g 1 -^ 2 ) 2 -(tany 1 -tany 2 ) 2 
4(^i + K 2 f - (tan - tan y/- 2 ) 2 ' 

tan = L ; /£T; , (/ = 1 , 2) . (26) 

which is the most well-known 2-soliton solution of the KP equation |[T8l . In the next 
section, we study more general soliton solutions which are expressed by the Wronskian 
form. 



3 Wronskian solutions, resonant-soliton and extreme 
elevations 

If the functions fi(X,Y,T) satisfy the relations 

dfi d 2 f; dfi d 2 fi 



dY dX 2 ' dT dX 3 



(i =!,-■•, AT), (27) 
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then the Wronskians 



(0) 



(N-l) 



AO) 
Jn 



JN-l) 

Jn 



with^ (n) = ^, (i = 1,..., N) (28) 



are solutions to the bilinear equation (fT9l |[T9l . For example, we can obtain the ordinary 
A-soliton solutions by setting /; as 

fi = e 9 »-i+e 9 », (i=l,-..,N), (29) 
8j = -kjX + k)Y -kjT + d®, (j= I,"- ,2A0, (30) 

where &y and are constants. 

Let us choose the following /; as satisfying (|2"7T ) 

A/ 

/i = X> e '=:/, f i =f {i - 1) , Ki<N<M,. (31) 



where 0,-, (/ = 1,2, • • • ,M) are given in the form (|30l and the ordering &i < k 2 < • • • < k M is 
assumed. Then T-function in (|28T) takes the form 



f(0) ... f(N 



-1) 



(32) 



y(JV-l) ... y(2AT-2) 

Using the Binet-Cauchy theorem, the T-function (l32l with (1311 ) can be expanded as a sum 
of exponential functions, 



£ A(/ b ...,/ M+ ) exp £ 0, 



A/ 



(33) 



i<ii<-<iM + <M v ;'=i 
where M = M + + M_, AT = AT+, (1 < M + < M— 1). A(/'i, . . . ,4?+) is the square of Van- 
dermonde's determinant, A(h, . . . ,iM+) = I1i</</<m + (^ ~~ ^'/) 2 lH21- This is called the 
fully resonant (M_,M + ) -soli ton solution. This solution has M_ solitons in asymptotics for 
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Figure 1 : Fully resonant soliton solutions 



y — > — oo, M + solitons in asymptotics for y — > °°, in the intermediate region these solitons 
interact resonantly. The T-function 031 is a general form of fully resonant-type soliton 
solutions. 



It is also possible to set as 

M 

fi = ^ aij e Q ' , for i = 1, . . . , N, and M > N , 

7=1 

using the phases dj given in the form d30l ) and a.N x M-matrix A( N ^ 
T-function (|28T ) with (1341 ) is expanded as a sum of exponential functions, 

t= £ £(h,...,i N ) 11 (*o~*i,)exp ( Y, ij ) . 

l<il<-</jv<Af i<j<l<N \j=l / 



(34) 



(a/,-). Then the 



(35) 



where £ (ii, . . . , /#) is the x A^-minor whose j-th columns (j = 1 , • • • , N) are given by the 
iy-th columns in the matrix Ar N<M \ = (aij), respectively, 



a i 



'A' 



§(ii,...,/jsr) 

This is a general form of solutions in which the line solitons interact resonantly and non- 
resonantly |[T3l . The T-function (f33l of the fully resonant line soliton solutions is a special 
case of the general T-function (T35T ). In this general form of line soliton solutions, 2-soliton 



solutions in the sense that the solutions have the same sets of two line solitons in both 
asymptoyics for y — ► ±oo are classified into three types, O-type, P-type and T-type. The 
A-matrices Aw^i corresponding these three types have the following row reduced echelon 
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Figure 2: The interaction pattern of the fully resonant (2, 3)-soliton {k\ = — 2,&2 = — 1, 
-1/100,^4 = 1/100,^ = 2) 



forms 



i(O) 
L (2,4) 



110 0^ 
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L (P) 

L (2,4) 



'l 



\ 



I 



L (2,4) 



/ 



V 



i o 
o 1 



+ 



110 

respectively, where '+,— ' shows the signs of the entries (also nonzero). T-type soliton is 
resonant two soliton solution, and O-type soliton is corresponding to two soliton solution 
obtained by the Hirota direct method lfT7l . Note that all jty-soliton solutions in the sense 
that the solutions have the same sets of N line solitons in both asymptoyics for y — ► ±oo 
are classified by using Young diagrams, i.e. we can construct the Any 2A0~ matr i ces from the 
Young diagrams |fT3l . 

Biondini and Kodama |[T2l and Kodama lfT3l considered the interaction patterns but did 
not the amplitudes of line solitons. Here, we consider the amplitudes generated by soliton 
interactions. First, we study fully resonant soliton solutions given by the T-function (|33T >. 
FigureQ]shows three examples of resonant soliton solutions. FigureQJa) shows the Y-shape 
resonant soliton solution (resonant (2, l)-soliton). The two incoming solitons interact reso- 
nantly and generate a high amplitude soliton. The amplitude of the high amplitude soliton 
generated by soliton resonance is four times as large as the amplitude of the incoming soli- 
ton. We see the amplitude of soliton is amplified by soliton resonance. Figures QJb) and 
[flc) show the resonant (2,2)-soliton (T-type 2 soliton). In general, T-type 2-soliton has a 
hole as seen in FigQJb). However, in FigQJc) we see no hole but H-shape soliton. This 
is a degenerate case of the T-type 2-soliton which occurs due to the closeness of two pa- 
rameters &2 and kj. The amplitudes of the intermediate interaction region in Fig{]Jb) are 
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Figure 3: The interaction pattern of the fully resonant (3, 3)-soliton {k\ = —5/2,^2 = 
-5/4,* 3 = -l/2,k 4 = l/2,k 5 = 3/2,fc 6 = 5/2) 

not so high, i.e. smaller than the double of the amplitude of the incoming soliton. On the 
other hand, in FigfTJc) we see the amplitude of the interaction region of the H-shape soli- 
ton is very near to four times the amplitude of the incoming soliton, i.e. higher than the 
sum of the amplitudes of two incoming solitons. In FigQTb) the interaction region makes a 
hole, and then the amplitude spreads out, then the amplitude becomes higher. In contrast, 
FigHJc) has no hole and the amplitude concenteates on small region. Physically, the ex- 
istence of the H-shape soliton is predicted by Miles by using symmetry argument of the 
Y-shape resonant soliton ifTOTl - The H-shape soliton is not a degenerate solution of 2-soliton 
solution. The resonant 2-solitons having a hole have been found in numerical simulations 
of the (2+l)-dimensional soliton equations (201I2TJ. It is natural that the H-shape solitons 
appear in physical situations such as ocean surface waves. The Y-shape soliton is a limit of 
the H-shape soliton. 

Figure |2] shows the interaction pattern of fully resonant (2, 3) -soli ton. In contrast to or- 
dinary nonresonant soliton solutions, the interaction pattern of fully resonant (2,3)-soliton 
shows complex nonstationary behavior. 

Figure [3] shows the interaction pattern of fully resonant (3, 3) -soliton. We see there 
exist 4 holes. In general, resonant (M_,M+) -soliton solutions given by the T-f unction (1331 ) 
generate web-like structure having (M_ — 1)(M+ — 1) holes lfT2l . 

The existence of web structure makes us expect the existence of nonstationary solutions 
having maximum amplitude for finite time. Thus, we have a following question: are there 
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Figure 4: Time evolution of the maximum amplitude of fully resonant (3, 3)-soliton. Exact 



Figure 5: The time eveolution of the partial resonant soli ton solution given by the T-f unction 
(1351) with the A-matrix A (2 5) i.e., 2 x 2-Wronskian, ft = I^. =1 e e *, f 2 = l| =4 e e * (fci = 
-2, k 2 = 0,k 3 = 2, k 4 = 2.01 , k 5 = 4) 

Figure 6: Plots of figure[5]in the plane y =constant at t = 1: (a) y=5.5, (b) y=4, (c) y=-2 

any exact solutions having characters of extreme waves? In Figf2] and Figj3j we see that 
the maximum amplitude of fully resonant multi-soliton solutions decrease at certain time 
by contraries. Figure @] show the time evolution of the maximum amplitude of Figj3j Fully 
resonant soliton solutions having web structure do not have characters of the extreme waves. 

Next, we consider the T-function (l33T > which is the general form of AMine soliton so- 
lutions of the KP equation. This solution includes rich variety of solutions and there are 
nonstationary solutions which are similar to extreme wave. Figure [5] shows an example. 
The solution of Figj5]was made by the 2 x 2-Wronskian solutions, the A-matrix is 



solution. 



/ 



\ 



1110 



A 



(2,5) 



(36) 



V 



1 1 



/ 



We see that the stem of Y-shape soliton and the other soliton interact near-resonantly. When 



1 -soliton cross with the stem of Y-shape soliton, the interaction region is amplified into high 



amplitude (See Fig J5j Figj6]and FigJT]). Thus, if there are several resonant solitons and each 



interaction parts of resonant solitons interact resonantly, huge waves are generated. Using 
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this principle, we can easily understand that there is a mechanism of generation of the 



extreme waves. 



In FigJU we show the interaction of inelastic 2-soliton given by the An 4) -matrix 



/ \ 
1110 



^(2,4) 



(37) 



\^ 1 1 J 

We see that the interaction pattern is a combination of two Y-shape resonant solitons at 
t < 0, i.e. there are two resonant stems. At t = 0, two resonant stems have a head-on 
collision, then a big amplitude wave is generated at t > (see Figl9land[T0~l). 

In FigfTTl we show the contour plots of solitons corresponding to Figf5]and FigJU Based 
on the asymptotic analysis which was discussed in lPT2l . each solitons are labeled by the in- 
dex [i, j] which represents a soliton made by two dominant exponentials e 9i and e G J . Taking 
the limit k.4 — > £3, £5 — > k 4 > in the soliton solution of the case of A( 2 ,5) i n FigJ51 we can 
recover the interaction pattern of the soliton solution of the case of A( 2 ,4) in Figj8] Because 
a new soliton resonance is made by the limit, the arrangement of solitons is changed. 

The mechanism of generation of extreme waves is as follows: (i) Several solitons are 
generated by some reasons, e.g. fast ferry traffic, tsunami, current and topography, (ii) Sev- 
eral generated solitons interact resonantly, (iii) When the interaction region of a resonant 
soliton interacts with one of the other resonant soliton, the amplitude will be 16 times of 
the amplitude of a soliton. (iv) After the interaction of the stem of resonant solitons, the 
amplitude decreases. 

This shows that the soliton resonance can be a mechanism of generation of extreme 
waves. The resonant interaction of several solitary waves may cause extreme waves higher 
than 4 times of amplitude of a solitary wave. This mechanism may happen in ocean; for 
instance, tsunami, freak wave and solitary killer wave may be amplified by this mechanism. 
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Figure 7: The change of maximum amplitude of figure [51 

Figure 8: The time eveolution of the inelastic 2- soli ton solution given by the T-function (1351) 
with the A-matrix A ( 2 4) ([37]) (k\ = -2,k 2 = 0,fc 3 = 2,fc 4 = 4) 

Figure 9: Plots of figure[8]in the plane y =constant at t = 1: (a) y=5.5, (b) y=3.5, (c) y=-0.5 

4 Conclusion 

The KP equation has rich mathematical structure and many interesting solutions. Finding 
physically interesting solutions of the KP equation is an important problem. In this paper, 
we have searched solutions having characters of extreme waves and found such solutions 
which are related resonant soliton solutions expressed by the Wronskian form. We have 
found there exists the solutions similar to properties of the extreme waves on the sea surface. 
The mechanism of generation of extreme waves is, (i) several solitons are generated by 
some reasons, (ii) several generated solitons interact resonantly, (iii) when the interaction 
region of a resonant soliton interacts with one of the other resonant soliton, the amplitude 
will be 16 times of the amplitude of a soliton, (iv) after the interaction of the stem of 
resonant solitons, the amplitude decreases. 

We acknowledge helpful discussions with G. Biondini, S. Chakravarty and Y. Kodama. 
KM also acknowledges support from the 21st Century COE program "Development of 
Dynamic Mathematics with High Functionality" at the Faculty of Mathematics, Kyushu 
University. 
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Figure 10: The change of maximum amplitude of figure [8l 

Figure 11: The contour lines of interaction patterns in Figj5]and Figj8] Each solitons are labeled 
by the index [i,j] which represents a soliton made by two dominant exponentials e 6 ' and e e L 
The upper graphs are corresponding to Figf5J the lower graphs are corresponding to Figj8j 
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